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ABSTRACT 

We have developed a one- dimensional thermochemical kinetics and diffusion 
model for Jupiter's atmosphere that accurately describes the transition from the 
thermochemical regime in the deep troposphere (where chemical equilibrium is 
established) to the quenched regime in the upper troposphere (where chemical 
equilibrium is disrupted). The model is used to calculate chemical abundances 
of tropospheric constituents and to identify important chemical pathways for 
CO-CH4 interconversion in hydrogen- dominated atmospheres. In particular, the 
observed mole fraction and chemical behavior of CO is used to indirectly constrain 
the Jovian water inventory. Our model can reproduce the observed tropospheric 
CO abundance provided that the water mole fraction lies in the range (0.25 — 
6.0) X 10~^ in Jupiter's deep troposphere, corresponding to an enrichment of 
0.3 to 7.3 times the protosolar abundance (assumed to be H2O/H2 = 9.61 x 
10~^). Our results suggest that Jupiter's oxygen enrichment is roughly similar 
to that for carbon, nitrogen, and other heavy elements, and we conclude that 
formation scenarios that require very large (> 8 times solar) enrichments in water 
can be ruled out. We also evaluate and refine the simple time-constant arguments 
currently used to predict the quenched CO abundance on Jupiter, other giant 
planets, and brown dwarfs. 

Subject headings: Jupiter, atmosphere; atmospheres, chemistry; abundances, at- 
mospheres; planetary formation 
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Introduction 



The water abundance in Jupiter's deep atmosphere provides important clues for solar 
system formation and evolu tion and reveals co nditions in the solar nebula at the time of 
giant-planet formation (e.g.. iLunine et al.ll2004( ). In addition, the planetary water inventory 



has important implications for the cloud structure, energy balance, thermal structure, and 
chemistry of the Jovian troposphere. Unfortunately, the deep water abundance is difficult 
to obtain by remote sensing methods because H2O is expected to cond ense near the ~5 
bar level in Jupiter's cold upper troposphere (e.g., see iTaylor et al.l 120041 ). and clouds and 
other opacity sources limit the depth to which infrared and other wavelength radiation can 
penetrate. The only in-situ measurement of the Jovian water abundance, by the Galileo 
Probe Mass Spectrometer (GPMS), pertains to a meteorologically anomalous "hot-spot" 
region characterized by low cloud opacity, low mixing ratios of condens able species^ high 
thermal emission, and a water abundance which i ncreased with depth (e.g . 



19981 : [Niemann et al.l Il998l : iRagent et all Il998l : ISromovskv et all Il998l : IWong et al.l 120041 ). 



Orton et al 



1996 



Thus, it is unclear whether the probe descended deep enough to sample the well-mixed 
water abundance below t he cloud base, and the Galileo probe value of H2O/H2 = 4.9 ± 1.6 x 
10~^ at the 19-bar level (IWong et al.l |2004|) is genera l ly considered to b e a lower lirnit for 



Jupiter's Q/H inven tory (e.g., iRoos-Serote et al.ll2004j : iTavlor et al.ll2004j : IWong et al.l 12004 



de Pater et al.ll2005l ). For this reason, chemical models must be used to determine the deep 
water abundance in Jupiter's atmosphere until further measurements become available, such 

'3 Juno mission or future deep-atmosphere probes (e.g.. 



Bolton et al. 


2006; 


Atreva 


2004) 



Several investigators have used clever methods to estimate the deep H2O abundance of 
the giant planets by considering the observed tropospheric abundance of CO and other trace 
constituents and by investigating how H2O che mistry and atmospheric transport can influ- 
ence the abundance of these trace species ( e .g., 



1985, 19881 



However, 



Lodders and Peg. 



Prinn anc 



1994 I2OO2I : iBezard et al. 



Barshav 



2002 



1977 



Feglev and Prinn 



Visscher and Fegleyll2005l ). 



Prinn and Barshayl ( 119771 ). and all subsequent modelers who used their kinetic 



assumptions have been shown to be incorrect ( 


Dean and Westmoreland 


1987; 


Yung et al. 


1988 




Griffith and Yelle 


1999 




Bezard et al. 


2002 




Cooper and Showman 


2006 


). In addi- 



tion, some of the transport time-s cale arguments in the earlier w orks have been shown 
to have inappropriate assumptions ( ISmithI Il998l ; iBezard et al.l |2002|) . These problems offset 
each oth er to an extent, such that models using ear li er assumptions yi e ld rea sonable re- 
sults (cf. IVisscher and Fegleyll2005t IBezard et al.ll2002l ). iFegley and PrinnI ( 1l988l ) found that 
H2O/H2 = (0.46 — 5.8) X 10~ ^ is co nsistent with CO kinetics and atmospheric mixing on 
Jupiter, whereas IBezard et al.l (120021 ) derived H2O/H2 = (0.34 — 15.3) x 10~^ using revised 
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(Paee et al. 


1989; 


Yune et al. 


1988; 


Smith 


1998) 



and improved CO observations. Further refinement of this estimate was not possible in the 
Bezard et al. analysis due to uncertainties in reaction kinetics, convective mixing rates, and 
the back-of-the-envelope time-scale arguments used to derive the H2O abundance. 

Here we attempt to improve the determination of the deep water abundance in Jupiter's 
atmosphere by taking advantage of recent updates in thermodynamic parameters and reac- 
tion rate coefficients and by using a numerical model to provide a more rigorous quantitative 
test of the simple kinetic vs. transport time-scale approach. With a numerical model, we im- 
plicitly solve the continuity equations for all tropospheric constituents, considering reaction 
kinetics and atmospheric transport. As a result, our model tracks the transition from the 
thermochemical regime in the deep troposphere (where chemical equilibrium is established), 
to a quenched regime in the upper troposphere (where chemical equilibrium is disrupted). In 
contrast to previous studies, we make no a priori selection of the reaction mechanism, nor of 
the rate-determining step for the chemical conversion of CO into CH4. Instead, we input a 
full suite of chemical kinetic reactions connecting the different relevant tropospheric species, 
and allow the dominant chemical pathways for the conversion of CO— )'CH4 to be identified 
from our model results. Furthermore, we make no assumptions about the mixing length 
scale but instead model atmospheric transport via diffusion for an assumed eddy diffusion 
coefficient profile. We explore the effects of the tropospheric water abundance on the chem- 
ical behavior of CO and other oxidized carbon gases and use the observed CO abundance to 
indirectly constrain the water inventory in Jupiter's deep atmosphere. 

The paper is organized as follows. We begin with a description of our chemical model and 
the CO chemical constraint in ^ We present our model results in ^ derive an estimate 
of Jupiter's deep water abundance using CO as an observational constraint, identify the 
dominant chemical pathways involved for CO— i'CIl4 conversion in Jupiter's troposphere, 
and discuss the chemical behavior of other oxygen-bearing carbon species. In §1] we discuss 
implications of our results for constraining Jupiter's total oxygen inventory and planetary 
formation scenarios, and conclude with a summary in ||5l 



2. Chemical Model 



2.1. Numerical Approach 



We use the Caltech/JPL KINETICS code flAUen et al.lll98lr ) to calculate the vertical 
distribution of atmospheric constituents in Jupiter's troposphere by solving the coupled one- 
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dimensional continuity equations as a function of time t and altitude z for each species: 

Pi - U (V. 



drii 

dt dz 



where is the number density (cm^'^), $j is the vertical flux (molecules cm~^ s^^), Pj is the 
chemical production rate (molecules cm~^ and Lj is the chemical loss rate (molecules 
cm~^ s~^) of species i. 

The continuity equations are solved using finite-difference techniques for 144 atmo- 
spheric levels, with a vertical resolution of at least twenty altitude levels per sc ale height. 



Jupit er's pressure-temperature profile is taken from Galileo entry probe data (ISeiff et al. 



19981 ) from 17 bar to the 22 bar, 427.7 K level (i.e., the maximum depth achieved by the 
probe before its destruction) and is exten ded to greater temperatu res and pressures along 
an adiabat using the method described by iFegley and PrinnI (119851 ). A zero flux boundary 
condition is maintained at the top (17.4 bar, 399 K) and bottom (12,650 bar, 2500 K) of 
the model, with thermochemical equilibrium being used to define the initial conditions. The 
relative abundances of the elements are thus defined a priori, and no mass enters or leaves 
the system. Although our model does not include rock-forming elements which may react 
with oxygen in Jupiter's deep atmosphere, we do c onsider the part ial re moval of oxygen into 
rock ( see §2.31 below) based upon the approach of iLodderd (120041 ) and IVisscher and Fegley 
(120051 ) so that our abundanc e calculations involving ox ygen are correct for higher altitudes 
(T < 2000 K; cf. Fig. 39 in iFegley and Lodderd Il994j ). Model calculations are performed 
until successive iterations differ by no more than 0.1 



70. 



2.2. Eddy Diffusion Coefficient 



In the context of our one-dimensional model, transport is assumed to occur by eddy 
diffusion, characterized by a vertical eddy diffusion coefficient Kzz- In the absence of a 
strong magnetic field or rapid rotation, the eddy diffusion coefficient can be estimate d from 
free-convection and mixing- length theories (jStonelll976l : iGierasch and Conrathlll985l ). using 
the scaling relationship 



K ^ 



wH 



H 



V I (2) 
where w is the characteristic vertical velo city over which cq r ivection operates (cm s~M i F is 



Jupiter's internal heat flux (5.44 W m ^, iHanel et al.lll98ll : jPearl and Conrathlll99ll ). k-Q is 



the Boltzmann constant, p is the atmospheric mass density (g cm ^), m is the atmospheric 
mean molecular mass (g), and Cp is the atmospheric specific heat at constant pressure (erg 



- 5 - 



K ^ g ^). The characteristic length scale over which mixing operates is assumed to be the 
atmospheric pressure scale height if, given by 



H 



mg 



where g is the gravitational acceleration. Eq. ([2]) yields Kz 
much of Jupiter's deep troposphere. 



(3) 

7 X 10*^ cm^ s~^ throughout 



However, strong magnetic fi elds or rapid ro t ation can affect conve ction and alter K. 



[e.g., iFlasar and GieraschI 119771 : IStevensonI Il979l : iFernando et al.lll99ll ). Rotation reduces 
Kzz, whereas magnetic fields ca n either increase or decrease K^z over the estimates from 
standard mixing- length theorv (jStevensonl Il979l). Consi dering rotation alone, we use the 



equations developed by iFlasar and Gierasch 



(119771 . Il978l ) to describe small-scale thermally 



driven convection in a rapidly rotating atmosphere. As outlined by these authors, the char- 
acteristic vertical velocity w for turbulent convection under these conditions is 



w 



w 



i-agSf' 



Ro 
sin A 



for Ro ^ I sin A I 
for Ro <^ I sin A I 



(4) 
(5) 



where a is the thermal expansion coefficient (= 1/T for an ideal gas), 5* is the static stability 
89 /dr (K cm~^), i.e., the radial gradient of the potential temperature 9 (and note that S < 
for free convection, as we are assuming here), A is the planetocentric latitude, and Ro is the 
Rossby number associated with the convective flow: 



Ro 



-agS) 



1/2 



2fi 



(6) 



where VL is the rotational angular velocity (s^^). Note that Eq. (jlj), which is valid for near- 
equatorial latitudes, is the same expression that is derived from mixing-length theory in 
non-rotating systems (see Eq. ([2]) above). At non-equatorial latitudes, Eq. ([5]) is appropriate, 
although it is convenient to replace S with an expression involving F, the internal heat flux 
— a measured quantity. From mixing-length theory. 



w 



pCpHS 

After equating Eq. (|5]) and Eq. ([7]) and solving for S and Ro in terms of F, we flnd that 



(7) 



K ^ 



pmCp{4D?' sin A)^/^ 



3/5 



for Ro <ti I sin A I 



(8) 
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Figure [T] shows the solution to Eq. ([8]) for non-equatorial latitudes over our model 
pressure range, although we should note that the condition Ro <^ \ sinA| begins to break 
down at high altitudes and low latitudes (e.g., Ro = \ sinA| for 20 bar at 3° latitude). Note 
that these K^z values are smaller than those predicted by standard mixing-length theory in 
a non- rotating atmosphere (see above); moreover, there is a very strong latitude dependence 
of Kzz at any particular pressure. As noted by Flasar and Gierasch (1977), such meridional 
variations in K^z could have some interesting consequences, including the production of 
meridional gradients in temperature and in the mixing ratios of quenched disequilibrium 
species like CO. However, these gradients may in turn drive meridional fluxes that act to 
homogenize the temperatures and abundances. It is therefore difficult to determine an 
appropriate Kzz profile to adopt for our models. 

The Galileo probe entered the Jovian atmosph ere in the Northern Equatorial Belt (NEB) 



near 6.5° N latitude, and the CO observations of (IBezard et al.l 120021 ) were also centered in 
the NEB near 9° N latitude. At 9° latitude, Eq. (jH]) implies i^^^ ~ 4 x 10^ cm^ s"^ at 2400 
bar, increasing to 1 x 10^ cm^ s~^ at 250 bar and 2.5 x 10*^ cm^ s^^ at 20 bar. Because of 
uncertainties in the appropriate values of K^z in Jupiter's deep troposphere, we adopt an 
altitude-independent Kzz value of 1 x 10^ cm^ s~^ for our nominal model, but we also explore 
the effects of a range of altitude-independent Kzz values from 4 x 10^ to 1 x 10^ cm^ s~^, as 
well as consider Kzz values that vary with altitude using Eq. (|8]). We assume that this range 
encompasses all plausible values for Kzz near the quench level. 



2.3. Atmospheric Composition 



Our model includes 108 hydrogen, carbon, nitrogen, oxygen, and phosphorus species 
that are subject to vertical transport and chemical production and loss. The measured 
relative abundances of these elements (see Table [1]) are used to define our initial conditions. 
The observed abundances of sulfur (as H2S) and the noble gases Ar, Kr, and Xe are shown 
in Tabled] for comparison. Atmospheric mixing ratios for He, CH4 an d NH3 are taken from 



1998 



GPMS and helium i nterferometer measurements: He/H2 = 0.1574 (Ivon Zahn et al. 
Niemann et aDll998h . and CH4/H2 = 2.37 x lO'^ and NH3/H2 = 6.64 x 10"^ JWong et al. 



20041 ). Our adopted phosphorus elemental abundance is based on a deep phosphine mixing 



ratio of PH3/H2 = 8 x 10 derived fr om Voyager IRIS (IKunde et al. 



19891), Galileo NIM S (llrwin et a 



and ground-based (IBezard et al. 



1998r) , Kuiper Airborne Observatory (IBjoraker et al.lll986[ ) 



2OO2I ) 5-/im observations of Jupiter 



1982 



Lellouch et a 



The Galileo entry probe measured a tropo spheric water abundance of H2O/H2 = 4.9 x 
10~^ near the 19-bar level (IWong et al.l |2004[ ). which corresponds to a subsolar 0/H ra- 
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tio. As was previously mentioned, it is unclear whether this value is representative of 
Jupiter's deep, well-mixed water inventory because the probe entered a meteorologically 
dry hot-spot region and because the H2O abu ndance increased nearly tenfold over the 
previous measurement at the 11-bar level (see IWong et al.l 120041 ) . Subsolar 0/H values 



are also consistent with 5-Atm observati ons from the Kuiper Airborne Observatory (KAO) 



(ILarson et al. 



1975 



Bioraker et al.lll986|). the Vovaaer Infrared Interferometer Spectrometer 



(IRIS ) ( Drossart and Encrenaj 1982 : Kunde et al. 1982 : Bioraker et al. 



1986 



1989h . the Galileo Near- Infrared Mapping Spectrom eter (NIMSl (llrwin et al.l 
1998h . and the Infrared Space Observatory (ISO) ( Encrenaz et al. 19961 : 



Lei 



1998 



ouch et al. 



Roos-Serote et al 



Roos-Serote et al. 



I999I . I2OO4J ). However, the interpretation of the observations is complicated by several fac- 
tors: (1) these 5-fim observations are weighted by thermal emission from the relatively dry 
and cloud-free hot-spot regions, (2) the observations are most sensitive to an altitude re- 
gion in which water is expected to condense so that the H2O mole fraction will be varying 
strongly with altitude, and (3) the results are highl y dependent on assumptio ns regarding 
cloud structure and extinction. As is discussed by iRoos-Serote et al.l (l2004j ). the Galileo 
NIMS observations are consistent over all with a deep 0/H ratio of 1-2 times solar (see the 
discussion below on the solar ratio) and lRoos-Serote et al.l ( 120041 ) state that "subsolar values 
of the 0/H ratio cannot be reconciled with the analyzed data," but this conclusion is also 
model dependent. Due to the difficulties in determining the 0/H ratio from remote-sensing 
observations and due to the ambiguities in interpreting whether the probe data are repre- 
sentative of the whole planet, we treat the oxygen elemental abundance as a free parameter 
in our model — indeed, the main free parameter that we are ultimately trying to constrain 
from our model-data comparisons — and we consider a wide range of H2O/H2 ratios in order 
to explore the effect of this parameter on tropospheric chemistry and the CO abundance. 

Table [1] also contains the elemental enrichment factors relative to the assumed protosolar 
abundances. Elemental abundances for the solar nebula (i.e., protosolar abundances) are 



taken from iLodders et al.l ( l2009l ). These abundances are slightly different from present-day 
solar photosph eric abundances because of gravitational settling of heavy elements in the Sun 
( jLoddersll2003l ) and represent the bulk elemental composition of the Sun and solar nebula. As 
in previous studies (jLodderj |2004J : IVisscher and Fegleyll2005l ). the protosolar H2O/H2 ratio 
is defined by taking the total oxygen abundance (SO/H2 = 1.21 x 10~^) and subtracting the 
portion that forms rock (Orock): 



OH2O — SO — Orock- 



(9) 



In a gas with a protosolar composition ( iLodders et al.ll2009l ). the formation of rock effectively 
removes ~20% of the total oxygen inventory. Throughout the following, we thus define the 
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water enrichment factor {Eu^o) as 



Eh,o = (H2O/H2) Jupiter 



solar ; 



(10) 



in which enrichments over "solar" refer to a protosolar abundance of H2O/H2 = 9.61 x 10~^. 
Using this definition, the Galileo probe measured a water abundance equivalent to 0.51x 
solar (-EH20 = 0.51) in Jupiter's troposphere. 

The relative abundance of H2 (which comprises 86% of the atmosphere) is calculated 
by difference using the observed mole fractions for He, CH4, NH3, PH3, and varying as- 
sumptions for the H2O mole fraction. From the adopted abundances of these major com- 
pounds, the initial chemical composition for all 108 species in the model is calculated 
along Jupiter's pressure-temperature profile by thermochemic al equilibrium. We use the 



NASA CEA (Chemical Equilibrium with Applications) code (iGordon and McBridd Il994[ ) 
fo r our equilibriu r n calculations, uti li zing thermod y namic parameters f r om th e compilations 
of ICurvich e^aP Jl989[ Il99ll . Il994l ): Ichasel f ll998h : ISurcat and Ruscid toO^ and other lit- 
erature sources. Figure [2] shows our equilibri um results for carbon chem istry in Jupiter's 
atmosphere, which are similar to the results of iFegley and Lodderj ( 1l994j . see their Fig. 17) 
with the exception o f H2CO and CH2: here we us e the H2CO enthalpy of formation from 



epti 

Gurvich et all ^9891 ll99lL Il994[ ) instead of lChasd (119981) based upon the recommendation 



of Ida Silva et al.l ( 20061 ) . and we calculate the mole fraction abundances of ^CH ? and ^CH2 
individually using updated thermodynamic parameters from lRuscic et al.l (120051 ). From this 
initial equilibrium condition, we run the kinetic-transport model to solve the continuity 
equations (described by Eq. [1]) until steady state is achieved. 



2.4. Chemical Reactions 



Approximately 1800 chemical reactions are included in our model. Photolysis reactions 
are omitted because the tropospheric pressure levels under consideration (> 17 bars) are too 
deep for ultraviolet photons to penetrate. Ou r reaction list is based upon previous chem - 



ical na od els of giant-planet atmospheres (e.g., iGladstone et al.l Il996l : iMoses et al.l ll995al Jbl. 



2000al Jbl. |2005[ ) and includes numerous r ate coefficient updates from combustion chemistry 



studies over the past two decades fe.g.. iBau 



Dean and BozzeUi 2000; Korobeinichev et al. 



chetal 



19921. 11994. 



200d iGardir^ 



2000 



2005 


Smith et al. 


1999 


Miller et al. 


2OO5I). 



All of the reactions in the model are reversed using the principle of microscopic re- 
versibility, 

k 



K = 



11) 
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where Keg is the equihbrium constant, kf is the rate coefficient for the forward reaction 
and kr is the rate coefficient for the reverse reaction. The the rmodynamic p a ramet e rs for 
calcu l ating K^^ are take n primarily from the com pilations of iGurvich et al.l ( 119891 . Il991 



19941 ). IChasd (Il998l ). and lBurcat and Ruscid (120051 ). We also adopt updated th ermodynamic 
data for sev eral radical species (notably OH, CH3, CH2OH and CH3O) from iRuscic et al. 



(120021 . 12OO5I ). Equation (fTTj) is used to determine kr for every reaction at each of the 144 
atmospheric levels using the appropriate temperature-dependent values for kf and K^q. The 
outcome of this procedure is a reaction list which consists of ~900 forward reactions and ~900 
complimentary reverse reactions, which can be used to accurately reproduce equilibrium 
abundances based upon the selected thermodynamic parameters. In other words, in the 
absence of disequilibrium effects such as photochemistry or atmospheric mixing, and given 
sufficient time to reach a steady state, our kinetic model results are indistinguishable from 
those given by Gibbs energy minimization or mass-action, mass-balance thermodynamic- 
equilibrium calculations. 

For three-body (termolecular) reactions, we assume that the forward reaction rate con- 
stant kf (cm^ s~^) obeys the expression 



ko 



ko[M] 

^00 



(12) 



where ko is the low-pressure three-body limiting value (cm^ s~^), k^o is the high-pressure 
limiting value (cm^ s~^), and [M] is the total atmospheric number density (cm~^). The 
exponent (3 in Eq. ( fT2|) is given by 



f3 



1 + 



log 



10 



fco[M] 



2 

J / 



(13) 



and we assume Fc ~ 0.6 (IDeMore et al.lll992l ) except where specified in lMoses et al.l ( 120051 ). 
Termolecular and unimolecular reactions in the model are also reversed using the method 
described above (Eq. ( ITT]) ). 



2.5. CO Chemical Constraint 



(1978 


), 


Larson et al. 


(1978 


), and 


" ■■ ' 

Noll et al. 



( 119881 . 11997I ). High-spectral-resolution observations that resolve the line shapes are needed 
to determine the origin of the CO, which can have either an internal source (i.e., from rapid 
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mixing from the deep troposphere) or an external source (e.g., from satelhte or meteoroidal 
debris). We are concerned solely with the internal source in this paper. As such, we rely 
on the observations that have most definitively resolved the CO vertical profile, those of 
Bezard et al.l (|2002[ ). Fro m ground-bas e d ^4 .7 /xm observations with a very high spectral 



resolution of 0.045 cm ^, iBezard et al.l ( 12002| ) determined that both internal and external 



sources contribute to CO on Jupiter and concluded that the derived tropospheric mole 
fraction of (1.0 ± 0.2) x 10~^ represents the contribution from the internal source. 

Our basic approach for const raining Jupiter's water inventory is similar to the method 

adop t ed by previous authors (e.g..lFeglev and Prinn ll988l : lFeglev and LodderJlQoi iLodders and Fegleyl 



I994J : iBezard et al.ll2002l : IVisscher and Fegley 



2OO5I ) : the observed abundance and chemical 



behavior of CO is used to estimate the H2O abundance in Jupiter's atmosphere. Deep in the 
troposphere, CO is produced from water and methane via the net thermochemical reaction 



CH4 + H2O = CO + 3H2 



(14) 



At high pressures and low temperatures, the reactants on the left-hand side of the reaction 
(CH4, H2O) are favored, whereas CO becomes more stable at low pressures and high tem- 
peratures. On Jupiter, the temperature factor is more important, and the CO abundance 
is expected to increase with depth, although methane remains the dominant carbon-bearing 
gas throughout the atmosphere. The equilibrium constant expression for reaction f|T4|) may 
be written as 



[CO][Hs 



(15) 



[CHjpaO]' 

where -ftjxi] equilibrium constant for reaction f|T^ and [i] is the concentration (cm^^ 

of each species i. Rearranging this expression. 



[CO] = i^^[CH4][H20]/[H2] 



(16) 



show s that [CO] oc [H9O] at e quilibrium for constant p ressure and temperature conditions 



Fegley and Primj Il988l : IVisscher and Fegleyl l2005l ). Thus, if thermochemical equilib- 



(e.g., 

rium holds, the abundance of any one of these species in reaction ( |T^ at a given altitude (i.e., 
at a given pressure and temperature) can be readily calculated from the observed abundances 
of the other three species at the same altitude in Jupiter's atmosphere. 

However , the observed CO mole fraction of Xqq = (1.0 ± 0.2) x 10~^, reported by 



Bezard et al.l (120021 ) for the 6-bar level in Jupiter's atmosphere, is many orders of magnitude 
higher than the CO abundance predicted under thermochemical equilibrium conditions for 
any plausible assumption about the deep water abundance on Jupi ter and provides clear ev- 
idence of disequilibrium processes at work in the troposphere (e.g 



Prinn and Barshav 



1977 



Barshay and Lewis 1978; Fegley and Prinnlll985 . 1988; Fegley and Lodderslll994 : Bezard et al. 
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20021 ). Carbon monoxide is not in equilibrium in Jupiter's upper troposphere, and therefore 
Eq. f|T6|) cannot be used to derive Jupiter's deep water inventory. 



Prinn and Barshayl ( 119771 ) demonstrated that the observed CO abundance hkely resuhs 



from mixing from deeper atmospheric levels where CO is more abundant. As parcels of 
gas rise in Jupiter's atmosphere, CO is destroyed by conversion into CH4 at a rate that 
falls off dramatically with decreasing temperatures. The observable amount of CO in the 
upper atmosphere thus depends upon the relative time scales of CO destruction kinetics 
(characterized by Uhern) and convective vertical mixing (characterized by tmix) (e.g., see 
Prinn and Barshayl 1 19 771 : iFegley and PrinnI 119851 ). Deep in the troposphere, tchem < imix 



and thermochemical equilibrium is achieved because kinetic reactions dominate over atmo- 
spheric mixing. At higher, colder altitudes, tchem > tmix, and disequihbrium prevails because 
vertical mixing operates faster than reaction kinetics can attain equilibrium. At some inter- 
mediate "quench" level defined by tchem = tmix, the CO mole fraction becomes quenched; 
above this level, the CO abundance rernains f ixed at the equilibrium mole fraction achieved 
at the quench level ( iPrinn and Barshayl 119771 ). As a result, the observed CO mole fraction 
in Jupiter's upper atm osphere still serves as useful chemical probe for conditions in Jupiter 's 
deep atmosphere (e.g.. lFegley and Prinnlll983l : lFegley and Lodderslll994j : lBezard et al.ll2002l ). 
despite the departure from equihbrium, provided that the kinetics of the CO — )■ CH4 con- 
version process is accurately known (to define tchem) and the rate of atmospheric mixing can 
be constrained (to define tmix)- 

Although we go beyond the back-of-the-envelope tchem-tmix approach with our kinetic- 
transport model, the underlying principle is the same. We determine the range of water 
enrichments that are consistent with the observed CO mole fraction for different plausible 
assumptions of the eddy diffusion coefficient and thereby indirectly constrain Jupiter's global 
water inventory. 



3. Model Results 

3.1. CO Profiles and Their Sensitivity to Model Free Parameters 

The CO vertical profiles resulting from our thermochemical kinetics and diffusion model 
are presented in Fig. [3] for various assumptions about the eddy diffusion coefficient and for a 
single assumption about the deep water abundance (Ix solar). Sho wn for compar i son is the 



observed CO mole fraction of Xqq = (1.0 ± 0.2) x 10~^ reported by lBezard et al.l (120021 ) for 
the 6-bar level on Jupiter; note that we have assumed that this value remains constant down 
to at least the few tens of bar region. The dashed gray line in Fig. |3] shows the predicted 
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equilibrium abundance of CO (calculated using the NASA CEA code), which increases with 
depth. In the deep atmosphere (P > 600 bars, T > 1100 K), the modeled CO abundance 
follows thermochemical equilibrium because the energy barriers for chemical reactions are 
easily overcome at these higher temperatures. As a result, the characteristic chemical time 
scale for interconversion between CO and CH4 is much shorter than the time scale for 
vertical mixing (tchem tmix), and equilibrium is readily achieved. However, as atmospheric 
parcels are transported to higher, colder altitudes, the chemical kinetic timescale becomes 
longer (i.e., the reactions become slower) relative to the mixing time scale, and the CO mole 
fraction can become quenched. In Fig. [3l divergence from the equilibrium profile occurs 
where vertical convective mixing begins to drive the CO abundance out of equilibrium and 
toward a constant quenched mole fraction. The transition region between the equilibrium 
and quenched regimes (where tchem ~ imix) is not abrupt but occurs over a range of altitudes 
approximately equal to one pressure scale height: for example, at K^z = 1 x 10® cm^ s~^, 
the modeled CO abundance diverges from equilibrium near the 640-bar level and resumes 
a constant quenched profile above the 200-bar level. Our numerical models confirm the 



pron 

analytic prediction of IPrinn and Barshayl (jl977l ) that the CO mole fraction remains fixed 



at altitudes above the quench region when the chemical reaction kinetics for CO — CH4 are 
extremely slow relative to the rate of vertical mixing {tchem ^ tmix)- 

Figure [3] illustrates that the "quenched" upper tropospheric CO mole fraction depends 
on the strength of convective mixing. For stronger convective mixing (i.e., larger K^z), CO 
is quenched deeper in the atmosphere where the CO mole fraction is larger. Conversely, 
weaker mixing results in lower quenched CO mole fractions. If the deep water enrichment is 
Ix solar. Fig. [3] shows that the CO observations are best reproduced for Kzz = 4 x 10® cm^ 
s~^, if we assume a constant Kzz profile with altitude. If the Kzz profile varies with altitude, 
as with the model represented by a black dashed line in Fig. [3l the results are controlled by 
the Kzz value near the quench level. 

The quenched CO mole fraction also depends on the assumed deep water abundance. 
Figure H] shows vertical abundance profiles for CO in Jupiter's atmosphere for our assumed 
nominal Kzz value of 1 x 10® cm^ s~^ over a range of water enrichments, including the 
abundance measured by the Galileo entry probe (0.51x solar). The gray lines indicate the 
equilibrium abundance of carbon monoxide for different water enrichments; divergence from 
the equilibrium profiles occurs where fast atmospheric mixing and slow reaction kinetics 
(relative to one another) quench the CO mole fraction. For any single Kzz profile, the 
pressure level at which quenching occurs is roughly independent of the H2O abundance, 
over the range of water enrichments considered here. Therefore, for all the models shown in 
Fig. m the quench level occurs at similar pressure and temperature conditions. As discussed 
above, for otherwise constant conditions (e.g., P, T, K, [CH4], [H2]; see Eq. [16]), the CO 
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abundance is linearly proportional to the H2O abundance in Jupiter's deep atmosphere (see 
also Fig. |5]). The quenched CO mole fraction thus increases as the assumed deep water 
abundance increases. For the assumption K^z = 1 x 10^ cm^ s~^, Fig. H] shows that the CO 
observations are best reproduced for a global water enrichment between 2-4 times solar. 

Table [2] and Fig. [5] summarize the sensitivity of the quenched upper tropospheric CO 
mole fraction to variations in the two main free parameters of our kinetic/transport model: 
the atmospheric convective mixing rate (characterized by K^z) and the tropospheric water 
inventory (characterized by Ey^^q). Taking the CO mole fraction as an observational con- 
straint, we can constrain the deep H2O abundance on Jupiter for a plausible range of eddy 
Kzz values. For example, our model solutions are represented by the solid lines in Fig. |5l 
The slope of each line is proportional to K^z and hence the vertical co nyective mixiiig rate . 



The observed upper tropospheric CO mole fraction of (1.0±0.2) x 10~^ (IBezard et al.ll2002l ). 
represented by the shaded area, constrains the range of Kzz and E-^^o values in our model 
solutions which are consistent with observations of CO in Jupiter's troposphere. 

For our nominal model with Kzz = 1 x 10^ cm^ s~^, our results show that an H2O 
abundance of 2.5 times the solar H2O/H2 ratio best reproduces the observed CO abundance. 
This enrichment corresponds to H2O/H2 = 2.4 x 10~'^ or XhsO = 2.1 x 10^^ in Jupiter's 
troposphere. If we further consider a range of plausible Kzz values from 4 x 10^ to 1 x 10^ cm^ 
s~^ based on the theories described in Section [27H we find that H2O/H2 ratios between 0.52 
and 5.20 times the solar ratio remain consistent with CO observations, including uncertainties 
in the observations themselves. 

The above discussion does not consider the potential errors due to uncertainties in 
reaction kinetics and thermodynamic parameters. The error bars in Fig. |5] show our at- 
tempt to quantify the effect of these uncertainties. As described below ( §3.2p . the CO — )• 
CH4 conversion is largely controlled by one rate-limiting step (reaction IR863p . for which 
the rate coe fficient was deriy e d usin g the rate coefficient of the reverse reaction flR862p . as 



reported by iJodkowski et al.l ( 119991 ) from ab initio calculations. However, no uncertainties 
were discussed for the calculated rate coefficient %862l Based upon a literature search of 
theory-data comparisons for other reactions, we estimate that AjRgel] is uncertain by a factor 
of ~3, which dominates over uncertainties in the thermodynamic parameters for the calcu- 
lation of A |Ft863 l - Including this estimated factor-of-three uncertainty in the reaction kinetics, 
our model results give a water enrichment of 0.3 to 7.3 times the solar abundance, corre- 
sponding to H2O/H2 = (0.29 - 7.0) X 10^3 g^n^j x^i^o = (0.25 - 6.0) x 10"^ Thus, the 



subso lar water abundance (0.51x solar) measured by the Galileo entry probe (IWong et al. 



20041 ) is plausibly consistent with the observed chemical behavior of carbon monoxide if rel- 



atively rapid vertical mixing rates (e.g., corresponding to Kzz ~ 1 x 10^ cm^ s ^) prevail 
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in Jupiter's troposphere. For comparison, iBezard et al.l ( 120021 ) derive a water abundance of 
H2O/H2 = (0.34 — 15.3) X 10""^ (0.4 to 15.9 times the solar abundance) using a timescale ap- 
proach rather than a full kinetic-transport model. Refinement of our estimate for the water 
enrichment may require updated rate-constant measurements for key reactions in the CO 
reduction mechanism and improved constraints on the Kzz profile in Jupiter's troposphere. 



3.2. Reaction Pathways for CO Destruction 



In the present study we did not select the rate-determining step a priori, but rather 
we input a full reaction list and let the code identify the predominant (fastest) chemical 
pathway for net CO— )-CH4 conversion based upon the relative rates of all reactions included 
in the model. As discussed above ( §2.4p . reaction rate coefficients have been obtained from 
experimental or theoretical data available in the literature or calculated from the rate coeffi- 
cient of the reverse reaction at each level in the model via Eq. f llip . Empirical rate-constant 
expressions for the important reactions involving CO reduction in our kinetic scheme are 
listed in Table [31 

Our model results indicate that carbon monoxide reduction to methane in Jupiter's 
atmosphere is dominated by the following series of reactions 



H + CO ^ HCO 
H2 + HCO ^ H2CO + H 
H + H2CO ^ CH3O 
H2 + CH3O ^ CH3OH + H 
H + CH3OH ^ CH3 + H2O 
H2 + CH3 ^ CH4 + H 

CO + 3 H2 ^ CH4 + H2O 



(R742) 
(R787) 
(R831) 
(R863) 
(R858) 
(R151) 

(net) 



where M refers to any third body (i.e., the reaction is termolecular) and the reaction num- 
bers are the assigned reaction numbers in our ki netic model (see Ta b le [3l). This kinetic 
scheme differs considerably from that proposed by IPrinn and Barshayl (119771 ) and adopted 
by subsequent authors, which assumes that the rate-limiting step for converting CO to CH4 
is the reaction 

H2 + H2CO ^ CH3 + OH. (R669) 



Prinn and Barshavl (119771) calculated a rat e constant for IR669I using the rate of the reaction 
CH3 + OH (jFenimordll969l : lBowmarull974l ). assuming that formaldehyde is a major product. 
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However, a number of experimental and theoretical kinetic investigations over the past three 



between methyl and hydroxyl radicals and that the rate estimated by I 


-"rinn and Barsha^ 


1 


(1977 


is incorrect 


(Bar- Nun and Podolak 


19851: Dean and Westmoreland 


19871 Yune et al 




IQSsI; 


De Avillez Pereira et al.lll997: 


Xia et al. 2001: Krasnoperov and Michael 2004: Baulch et al. 


2005 


Jasper et al. 


2007), as has been noted in other studies of CO quenching kinetics 


(Yune et al. 1988: ( 


Griffith and Yelk 


1999: 


Bezard et al. 2002: Cooper and Showman 2006) 





Updated second-order temperature-dependent rate co efficients for the reverse r e action OH 



-|- CH.s — )• H2CO -|- H2 (R668) are now available from iDean and Westmoreland! (119871 ) and 



De Avillez Pereira et al.l (Il997( ). and show t hat the rate coeffici e nt for IR669I is ~3 orders of 



magnitude smaller than that estimated by iPrinn and Barshayl (119771 ) at quench-level tem- 
peratures (~1000 K) on Jupiter. Furthermore, we note that IR6691 is included in our model 



ics literature (e.g., I 


^enimore 


1969: 


Dean and Westmoreland 


1987: 


De Avillez Pereira et al. 


1997: 


Xia et al. 


2001 


: Krasnoperov and Michael 


2004 


Baulch et al. 


2005: 


Jasper et al. 


2007 





because faster, alternative reaction pathways exist for C0^CH4 conversion in Jupiter's tro- 
posphere. 



Our rea ction schera e (in p articular reaction s IR831I through IR151D i s similar to that 
pro posed bv lYung et al l Jl988l ) and adopted by ICriffith and Yellel Jl999l ) for Gliese 229B 
and bezard et all (|2002[ ) for Jupiter. However, we find that H2 + CH3O CH3OH + H 
(reaction IR863P is th e slowest reaction i n this dominant scheme and, as such, is the rate- 
limiting step, whereas lYung et al.l (Il988[ ). ICriffith and Yellel (119991 ). and lBezard et al.l (|2002[ ) 
assumed H + H2CO + M CH3O + M ( IR83ip was the rate-hmiting reaction for CO 
destruction. We also note that the activation energy for IR863lat high temperatu r es ( ~90 k J 



mol~^ for 1600 to 360 K) i s similar to that measured by IBar-Nun and Shavivl (119751 ) and 
Bar-Nun and PodolakI (119851 ) for CO reduction in high-temperature shocks. However, we are 
unable to confirm whether CO— >-CH4 conversion in the shock-tube exper iments involves th e 
same reaction pathways (IR7421 through [R151I) as we find for Jupiter (cf. lYung et al.lll988l ). 



There are no experimental data available for IR8631 so we have determined its rate from 
the rate coefficient of the reverse reaction. 



H + CH3OH ^ CH3O + H2, 



(R862) 



as calculated by IJodkowski et al.l (119991 ) using transition-state theory. A dopting hiKffa \ as fc r 
in equa tion (flTI) alori g with thermodyn amic paramet e rs for CH3OH from lChen et al.l (Il977l ). 
H from IChasd (Il998l ). and CH3O from iRuscic et al.l (120051 ) . we calculated the reaction rate 
coefficient for 1R8631 at each temperature level in our model and determine an empirical fit 
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of the form 



%863]=l.r7xl0-22T3°V-3055/T), 



(17) 



We note that IJodkowski et ahl (Il999l ) also calculated A |r863 | from the rate coefficient of the 



Jodkowski et al. 



revers e reaction, and estimate ^ggs] = 2.10 x iQ--^^x e~ ' . However 
(119991 ) use equilibrium constants derived theoretically from molecular parameters, wherea s 
we have taken advantage of the thermodynamic updates provided by Ruscic et al.l ( 120051 ) . 



Our value for A ]r863 I therefore differs from that of IJodkowski et al.l (Il999l ). but by an amount 



that is much smaller than our overall factor-of-three estimated uncertainty in A |Rgg3 | over the 
range of temperatures considered in our model. Our results regarding the quench point, 
and its constraint on the water enrichment, are particularly sensitive to the rate coefficient 
adopted for IR863t and the factor-of-three uncertainty in ^[^353] has been included in the 
overall uncertainties in the deep water abundance we derive for Jupiter (see §3.1 and the 
error bars in Fig. [5]). However, there are several alternative reactions, including IRSSlt which 
may dominate if this rate is in serious error, in which case we would expect qualitatively 
similar results for CO quench chemistry in Jupiter's troposphere (see note added in proof). 



3.3. Validity of the Time-scale Approach 



Having identified the dominant chemical mechanism for CO — t- CH4 conversion in the Jo- 
vian troposphere, we can now test the validity of the time-scale approach previously used for 



op 

estiraating Jupiter's deep water inventory (e.g., iPrinn and Barshaylll977t iFegley and Prinn 



19881 : iBezard et al.ll2002[ ). The chemical lifetime for CO is given by the expression 

[CO] [CO] 



chem 



(CO) 



-d[CO]/dt %863l[H2][CH30]' ^ ' 

assuming reaction IR863I is the rate-limiting reaction. The vertical mixing time scale is given 
by 

-li- 

tmix "fT (-^^) 
J^zz 

where Kzz is the eddy diffusion coefficient and L is the characteristic length scale over which 
the mixing operates. The atmospher ic pressure sc ale height is traditionally used for L for 
these types of calculations. However, ISmithI ( 1l998l ) demonstrated theoretically that L ^ H 
is not appropriate and may lead to over-estim ates of t he mi xing length and the mixing time 
scale. Using the procedure recommended by ISmithI (1l998[ l. we obtain L ~ 0.12H for CO 
quenching kinetics on Jupiter. 



We have calculated the CO abundance at the quench level (i.e., where tchem = t. 
over a range of water enrichments (0.51x to 8x) and K^z values (4 x 10"^ to 1 x 10^ cm 



mix J 

9 
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s~^) for compariso n with our kinetic model results. For Kzz = 1 x 10^ cm^ s~^ and Xqq = 
(1.0 ±0.2) X 10~^ ( iBezard et al.ll2002[ ). our time-scale approach yields a water enrichment of 
2.9 ± 0.6 times solar, whereas our nominal kinetic-transport model yields -EH20 = 2.5 ± 0.5 
(not including uncertainties in K^z and /^Rges])- Given the overall additional uncertainties in 
atmospheric mixing rates and reaction kinetics, we conclude that the back-of-the-envelope 
time-scale approach gives a reasonably accurate estimate of the H2O ab u ndanc e on Jupiter, 
provided that the vertical mixing length scale L advocated by ISmithI (Il998l ) is used and 
that the appropriate rate-limiting reaction and rate coefficient are considered. Regarding 
the latter point, we reemphasize that IR863] is the rate-limiting reaction for CO destruction 
in our model, in contrast to the assumptions of previous investigators. 



3.4. Chemistry of Other C-O Gases 

Carbon monoxide is not the only species that will undergo quenching in Jupiter's atmo- 
sphere. In principle, any atmospheric constituent subject to vertical transport and reaction 
chemistry will quench if the characteristic time scale for convective mixing becomes shorter 
than the characteristic time scale for kinetic destruction. Here we examine the chemical 
behavior of other oxidized carbon gases in Jupiter's troposphere. Figure |6] shows the vertical 
abundance profiles for C-O species for a model in which Kzz = 1 x 10^ cm^ s~^ and -EH20 = 1- 
The dotted gray lines indicate abundances predicted by thermochemical equilibrium using 
the NASA CEA code; divergence from equilibrium is evident for CO, CO2, H2CO and CH3O 
and illustrates where rapid vertical mixing and slow reaction kinetics (relative to one another) 
drives each species toward a constant quenched mole-fraction profile. Again, CH4 and H2O 
remain the dominant carbon- and oxygen-bearing gases, respectively, throughout Jupiter's 
troposphere. The CO abundance begins to diverge from equilibrium near the 640-bar level 
and assumes a constant quenched profile at altitudes above the 200-bar level. 



As pointed out by iPrinn and Barshayl ( 119771 ). reactions among oxidized (CO, CO 



^2, 



etc.) or reduced (CH4,CH3, etc.) carbon-bearing gases are expected to be much faster than 
reactions between these two families. As a result, the chemistry of many oxidized carbon 
gases is strongly tied to the chemical behavior of CO, the most abundant C-O compound 
in Jupiter's atmosphere. For example, CO quenching (via the rate-limiting step IR863|) 
immediately affects the vertical abundance profile of CO2. At altitudes above the level where 
CO begins to quench (~640 bar), CO2 remains in equilibrium with CO via the reactions 

OH + CO ^ CO2 + H (R708) 
H + CO2 ^ CO + OH . (R709) 

The CO2 abundance thus slightly increases with altitude until C0<=^C02 conversion (via 
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R708<=^R709) itself quenches near the 500-bar level, whereupon CO2 assumes a constant 
ve rtical profile to higher altitudes. Our upper-tropospheric result for CO2 is similar to that 



of iLellouch et al.l ( 120021 ). who used a time-scale approach to derive a quenched CO2 mole 
fraction abundance of 3 x 10^^^ for a solar 0/H ratio in Jupiter's troposphere. 

Although methanol is one of the products of the rate-limiting reaction (]R863P in our 
predominant kinetic scheme for CO— )'CH4 destruction, the vertical abundance profile for 
CH3OH shows no dependence upon the chemical behavior of CO. This result can be explained 
by comparing the most important reactions for CH3OH production at the 400-bar level (for 
example) in our nominal model: 

H2O + CH3 ^ CH3OH + H 93.5% (R859) 

H2 + CH2OH ^ CH3OH + H 6.3% (R861) 

H2 + CH3O ^ CH3OH + H 0.3% (R863) 

Methanol production from CH3 and H2O via R859 far exceeds production via R863, so the 
chemical behavior of CH3OH is largely decoupled from the 00—7-0134 kinetic scheme (i.e., the 
carbon in CH3OH comes predominantly from CH4— )-CH3 rather than CO, and the oxygen 
comes from H2O). The CH3OH abundance follows an equilibrium profile until production via 
R859 quenches near the 200-bar level, at which point vertical mixing drives CH3OH toward 
a constant mole-fraction profile. 

The vertical abundance profile for H2CO is more complex than that for the other 0-0 
gases. At altitudes above the level where CO begins to quench, HCO remains in equilibrium 
with CO via IR742t =^R743 (see Table [3]). In turn, the H2CO abundance departs from its 
predicted equilibrium abundance profile and instead remains in approximate equilibrium 
with CO because formaldehyde production via the reactions (at the 400-bar level): 

H2 + HCO ^ H2CO + H 82.5% (R787) 

CH3O ^ H2CO + H 17.5% (R830) 

is balanced by destruction via the reactions: 

H2CO + H ^ HCO + H2 81.1% (R786) 

H + H2CO ^ CH3O 18.9% (R831) 

This balance continues until conversion via R830<=^R831 begins to quench near the 500 
bar level and can no longer offset production and loss via HCO (R787^R786), whereupon 
vertical atmospheric mixing drives H2CO toward a constant quenched profile near the 100- 
bar level (a difference in altitude of roughly one pressure scale height). As a result (see 
Fig. [6]), the H2CO vertical abundance profile shows changes in response to two separate 
quenching reactions (R863 and R830) in Jupiter's troposphere. 
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Discussion 



4.1. Implications for Jupiter's Total Oxygen Inventory 

Water vapor is the dominant oxygen-bearing gas throughout Jupiter's atmosphere and 
is much more abundant than other oxygen-bearing gases (CO, OH, etc.). For this reason, 
the H2O abundance in the troposphere is expected to be representative of the majority 
of Jupiter's total oxygen inventory. However, some oxygen (^20 %) is removed from the 



gas phase by oxide formation (rock) in the deep atmosphere (e.g.. iFegley and Prinrull988 
Lodderd l2004l : IVisscher and Fegleyl 120051 ) . This removal must be considered before evaluat- 
ing the bulk planetary oxygen abundance and, in turn, the heavy-element composition of 
planetesimals during Jupiter's formation. 

The fraction of oxygen removed by rock depends upon the abundance of all rock-forming 
elements (Mg, Si, Ca, Al, Na, K, Ti) relati ve to the total oxygen abundance. FoUow- 
ing the method of IVisscher and Fegleyl (120051 ) . and using updated solar abundances from 
Lodders et al.l (120091 ) . the relative abundances of water vapor, total oxygen, and the rock- 
forming elements can be written as 



Eh,o = 1.261Eso - 0.261E,ock 



(20) 



where Ei represents the enrichment (over solar ratios) for each component i. This expression 
serves as a general mass-balance constraint for the relative abundances of water, oxygen, and 
rock over a range of heavy element enrichments in Jupiter's interior (IVisscher and Fegley 
20051 ). To derive this expression, it was assumed that all of the rock- forming elements are 
equally enriched, and -Erock = 2.74 =b 0.65 was adopted based upon the "deep" tropospheric 
abundance of sulfur (as H ^S jWong et al.ll2004l ). which behaves as a rock-forming element in 
meteorites (lLoddersll2004l ). Using E^ock = 2.74 ±0.65 in equation (l20l) along with the Galileo 
entry probe H2O abundance of E^i^o = 0.51 ± 0.17 (see Table [1]) yields a total oxygen 
abundance (characterized by -Eso) of 0.97 ± 0.40 times the solar abundance of SO/H2 = 
1.212 X 10^3_ 

Using our nominal model result of En^o = 2.5 ± 0.5, equation (!20l) yields a total oxygen 
enrichment of Ej^o = 2.5 ± 0.8. Further considering a 3x uncertainty in reaction kinetics 
and a range of plausible Kzz values from 4 x 10^ to 1 x 10^ cm^ s~^, our water constraint 
(0.3 - 7.3x solar) gives a total oxygen inventory of 0.7 to 6.5 times the solar SO/H2 ratio 
in Jupiter's interior. This value represents the bulk oxygen inventory of Jupiter's interior 
consistent with CO chemistry, and includes oxygen as water plus oxygen bound in rock. 
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4.2. Implications for Planetary Formation 



All viable giant-planet formation models must consider how heavy elements become 
entrained in the planet during its formation and evolution. An important constraint for 
such models is therefore the observed atmospheric abundances of gases such as CH4, NH3, 
H2S, PH3, and H2O, which are taken to represent the planetary elemental inventories of C, N, 
S, P, and the majority of planetary oxygen, respectively. As illustrated in Fig. [TJ observations 
of Jupiter's atmosphere sho w that C, N, S, P, Ar, Kr, and Xe are enhanced relat ive to solar 
element-to-hydrogen ratios (IMahaffy et al.ll200Cll : IWong et al.ll2004l : lLoddersll2004j ) by factors 
of 2-4. The enrichment in heavy elements is gen erally believed to be consistent with the 
core-accretion model for giant planet formation (IMizund 119801 ). in which a rock or rock- 
ice core initially forms and continues to grow thr ough the accretion of solid planetesimals 
until i t is massive enough to capture nebular gas (iBodenheimer and PoUackl Il986l : iLissauer 
19871 : iPoUack et al.l Il996l ). In this scenario, the observed heavy element enrichments on 
Jupiter arise from degassing of the initial core material and the continued acc retion of solid 
plane tesimals, which will most likely vaporize before reaching the core (e.g., iPoUack et al. 



19861 ). What remains unclear is the source and composition of the planetesimals which 
provided the enrichment, and several scenarios have been proposed to explain the observed 
heavy-element abundances. 

One distinguishing characteristic of Jovian formation scenarios is the predicted water in- 
ventory in Jupiter's deep atmosphere. For example, trapping of heavy elements by hypotheti- 



cal so lar-composition icy planetesimals ( Owen et al.lll999l : lAtreya et al.ll2003l : lOwen and Encrenaz 



20061 ) would be expected to give an enrichment in oxygen (as H2O) around 3±1 times so- 



clathrate hydrates n ear the snow line (e.g. 



Lunine and Stevenson 


1985; 


Gautier et al. 


2001a 



Hersant et al.ll2004r) would vield higher water abundances, with predicted enri chments rang- 
ing fro m ~6 ( Mousis et al.l2009) to ~ 8 ( Alibert et al. 2005 ) to as high as >17 ( Gautier et al. 
2001al lb[) or ~19 Jflersant et al.lbooi ) times the solar H-^ O/H^ ratioof 9 61 x 10""^. Accretion 



of carbon-rich planetesimals behind a nebular "tar line" (ILoddersI 120041 ) would give subsolar 
water abundances similar to that observed by the Galileo probe (0.51x solar). 

Our indirect constraint of a deep water abundance 0.3-7.3 times the solar H2O/H2 
ratio from our kinetic-transport model is plausibly consistent with each of these formation 
mechanisms but precludes clathrate- hydrate sce narios that would require large (> 8x) water 



enrichments in Jup iter's deep atmosphere (e.g.. iGautier et al.l 12001 
Alibert et al1l2005[ l. 



b; Hersant et al. 2004 
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5. Summary 



We have developed a comprehensive thermochemical kinetics and diffusion model for 
Jupiter which correctly transitions between equilibrium chemistry in the deep troposphere 
and quenched/disequilibrium chemistry in the upper troposphere. We use this numerical 
model to compute the vertical abundance profiles for all carbon- and oxygen-bearing atmo- 
spheric constituents and to explore the chemical behavior of CO and other oxidized carbon 
species in Jupiter's deep atmosphere. We find that ca r bon m onoxide is reduced to CH4 via 
a mechanism similar to that proposed by lYung et al.l (119881 ): however, our model indicates 
that the rate-limiting reaction for CO reduction in Jupiter's atmosphere is H2 + CH3O — )■ 
CH3OH + H rather than Yung et al.'s proposed re action H + CH3O + M -> CH3OH + M. 
We also confirm the original analytic prediction of iPrinn and Barshayl (Il977l ) that the mole 
fraction of CO will "quench" and remain constant with altitude when kinetic reaction rates 
can no longer compete with atmospheric mixing. This quenching occurs at the ~400 bar 
(1000 K) level in our nominal model. Carbon monoxide is not the only species to quench; 
virtually all atmospheric constituents will quench at some point where temperatures become 
low enough to inhibit the kinetics. 

Our kinetic-transport model quantitatively confirms the convenient, back-of-the-envelope 
time-scale approach currently used to explore quenched disequilibrium chemistry on giant 
planets and brown dwarfs (e.g.. 



1985 



19881 : iLodders and Fegley 



rinn and Barshav 



1994 



1977 



Lewis and Feglev 



2OO2I : [Griffith and YellelE999l : 



1984 



Bezard et al. 



Feglev and Prinn 



I2OO2I : IVisscher and Fegley 



2OO5I ). We find that the the time-scale approach is valid for estimating Jupiter's water in- 
ventory, provided that the correct rate-limiting reaction is considered (which we find to be 
reaction R863, H2 -|- CH3O — )■ CH3OH -f- H) and proy ided that the mixing length L is 



determined via the procedure advocated by ISmithI ( 119981 ) 



Using the CO abundance reported by iBezard et al.l (120021 ) as our observational con- 
straint, our model-data comparisons indirectly constrain the Jovian deep water abundance 
to lie in the range 0.3-7.3 times the solar H2O/H2 ratio of 9.61 x 10~^. Our results sug- 
gest that the enrichment for oxygen (as H2O) is similar, to within uncertainties, as that 
for carbon, nitrogen, and other heavy elements — giant-planet formation scenarios that re- 
quire very large (> 8x) enrichments in the water abundance (such as some clathrate-hydrate 
formation scenarios) are prec luded. The subso lar water abundance (0.51x solar) measured 
by the Galileo entry probe (IWong et al.l |200J) remains plausibly consistent with the ob- 



served tropospheric abundance of carbon monoxide if relatively rapid vertical mixing (e.g., 
Kzz ^ 1 X 10^ cm^ s~^) prevails in Jupiter's deep troposphere. 

We will not be able to narrow our estimated range of Jovian deep water enrichments 
without experimental confirmation of the rate coefficient for the reaction of methoxy with 
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H2 (CH3O + H2 — 7- CH3OH + H) at high pressures, for temperatures near 1000 K, and 
our results are subject to revision as updated kinetics data become available. Perhaps more 
importantly, we need a better understanding of appropriate diffusion coefficients that can be 
used to represent convective processes under tropospheric conditions on Jupit er. The Juno 
mission may provide microwave data ( Janssen et al. 2005 : Bolton et al.l l20o"6l ) that can be 



used to test our model prediction regarding the Jovian deep water abundance. 

We point out that oxygen species are not the only constituents to quench in our model; 
the quenching of nitrogen species like N2 and HCN is interesting in its own right and will 
be the subject of a future investigation. Our kinetics-transport model can easily be applied 
to other giant planets and brown dwarfs. Of particular interest is (1) constraining the deep 
water abundance on the other giant planets in our own solar system from current or future 
tropospheric CO observations (which must be able to constrain the vertical profile to separate 
the contributions arising from possible internal and external sources), as no current plans 
to send multiprobe missions to these planets are on schedule for the near future, and (2) 
predicting the vertical variation of observable species in brown dwarfs and extrasolar giant 
planets, as these atmospheres are unlikely to be in complete thermochemical equilibrium. 
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Note added in proof- In our kinetic mechanism for CO — ?■ CH4 conversion in Jupiter's 
troposphere, we adopted a rate coefficient for the re action H + CH3OH — t- CH3 + H2O 
(R858) based upon the work of iHidaka et al.l (Il989l ). However, recent literature studies 
suggest that this rate coefficient is inappropriate and may lead to over-estimates for the rate 
of (R858). The three main reaction pathways for H + CH3OH are: 



H + CH3OH ^ CH2OH + H2 
H + CH3OH ^ CH3O + H2 
H + CH3OH ^ CH3 + H2O 



(R860) 
(R862) 
(R858) 



Literature values for the re l ative rates of each pathway are contradictory (ILendvay et al. 



19971 : iJodkowski et al.lll999l : iBaulch et al.ll2005l ). Laboratory and theoretical investi gations 



indic a te that (R860) is the do r ainan t pathway with a rate coefficient that i s 4 times (ITsang 
198?! : iNorton and Drveil Il990l . Il989h to ~ 30 times JLendvav et al.l Il997l : Ijodkowski et al 
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19991 : ICarvalho et al.l |2008[ ) greater than that of (R862). These st udies also suggest that 



(R860) dominates over (R858) by a factor of ^10-40 at 1000 K fe.g..lAronowitz et al.lll977 



Hoyermann and Wagner! 1 19811 : ISpindler and Wagner Il982l : [Norton and Dryerlll990l ). except 
for lLendvav et all Jl997l ). who predict fcRgeo/^Rsss > 5000 at 1000 K. Althoug h (R858) is the 
most exothermic of the three reactions listed above (see al solYung et al.lll988l). this pathway 



appears to be inhibited by a large activation energy (e.g., iLendvay et al. 



19971 ) 



The rate-limiting step (RLS) for CO — )■ CH4 (and the chemical behavior of methanol) 
in Jupiter's troposphere is sensitive to both the overall rate of H + CH3OH and the relative 
contribution of each reaction pathway. We tested the sensitivity of our model results to 
variations in metha nol kinetics by adopting a rate coefficient for (R858) using fcRseo/^Rsss = 



40 at 1000 K 

-1 



mol 



le.g., 



[e.g., 



Norton and Drverlll990( ). along with an activation energy of ~ 103 kJ 



Lendvay et al.lll997l ). Depending upon the overall reaction rate and the relative 



rates of each pathway, we identify four possible rate-limiting reactions from our model results: 



H2 



CH,0 ^ CH,OH + H 



H + CH3OH ^ CH3 + H2O 



M 



H + H2CO A CH3O 



M 



CH3OH A CH3 + OH 



(R863) 
(R858) 
(R831) 
(R667) 



We note that (R863), (R858), and (R831) are each present in our mechanism for CO de- 
struction in Jupiter's atmosphere, and that the (R831) pa thway has been adopted as the rate- 
limit i ng step in previous s tudies of CO quenching kinetics ( Yung et al.lll988 ; Griffith and Yelle 
I999I : iBezard et al.l 120021 ). The reaction pathway (R667) will become the rate-limiting step 
only if the rate of (R858) is effectively negligible (R858 remains the RLS even at relatively 
high A;r86o/^R858 ratios, e.g.. iLendvav et al.l 119971 ). Considering a CO mole fraction abun- 
dance of 1.0 ±0.2 X 10-9 teezard et al.ll2002 h and a range of K^z values from 4 x 10^ to 1 x 10^ 
cm^ s"^, our model results using different H + CH3OH rate-coefficie nt data are summarized 
in Table |5l For our preferred model, we take fcRseo and /cr862 from IJodkowski et al.l (119991 ) 
and adopt A;r858 ~ 9.41 x io-9e(-^24oo/T) ^^^-3 g-i rYable 3), assuming fcRgeo/^Rsss = 40 at 
1000 K. In this approach, (R863) remains the rate-limiting step and the observed CO abun- 
dance (1.0 ± 0.2 ppb) yields a Jovian water abundance of 0.4-3.4 times the solar H2O/H2 
ratio (9.61 x 10^^), not including uncertainties in reaction kinetics. 

In each case, the water abundance derived from CO quench chemistry is consistent 
with our earlier results (which include an estimated factor-of-three uncertainty in reaction 
kinetics) of 0.3-7.3 times the solar H2O/H2 ratio. The different rate-coefficient data give 
roughly similar results because other reactions such as (R831) or (R858) may dominate if 
(R863) is not the rate-limiting step. Each of these reactions (R831, R858, R863) quench in 
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the same vicinity (within a fraction of a scale height) in Jupiter's troposphere. Our overall 
conclusions regarding the Jovian water inventory, and the reactions given in our kinetic 
scheme, thus remain unchanged. We emphasize the need for laboratory measurements of all 
the pathways for H + CH3OH in order to definitively identify the ratc-hmiting step for CO 
— )■ CH4 in Jupiter's troposphere and to refine estimates of Jupiter's deep water abundance. 
We gratefully acknowledge Greg Smith for alerting us to the uncertainties associated with 
our adopted rate for the CH3 + H2O reaction pathway. 
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Table 1: Gas Abundances in Jupiter's Troposphere 



Gas i 



Jupiter i/H2 



Reference 



Protosolar i/H2^ 



Enrichment Factor 



H2 
He 
CH4 
NH3 
H2O 
PH3 
H2S 
Ar 
Kr 
Xe 



0.1574 ±0.0036 
(2.37 ±0.57) X 10^3 
(6.64 ± 2.54) X 10"'' 
(4.9 ± 1.6) X 10-^ 
(0.8 ±0.5) X 10-6 
(8.9 ±2.1) X 10"^ 
(1.82 ±0.36) X 10-5 
(9.30 ± 1.70) X 10-9 
(8.90 ± 1.70) X 10-1° 



b, c 
d 
d 
d 
e 
d 
f 
f 
f 



= 1 
0.1938 
5.55 X 10-"^ 
1.64 X 10""^ 
9.61 X 10"^ 
6.41 X 10-^ 
3.25 X 10-5 
7.16 X 10-6 
4.31 X 10-9 
4.21 X 10-1° 



= 1 
0.81 ±0.02 
4.27 ± 1.03 
4.05 ±1.55 
0.51 ±0.17 
1.25 ±0.78 
2.74 ± 0.65 
2.54 ±0.50 
2.16 ±0.39 
2.11 ±0.40 



No te. References: (a) Pro toso lar abundances fro m lLodders et al.l ( 120091 ) fb MNiemann et al.l (119981) 
(c) Ivon Zahn et al.l (119981 ). (d) IWong et al.l (120041 ) . (e) 5-/im spectroscopy, see text (f) lMahaffy et al 
(|2OO0l ). 
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Table 2: Quenched CO Mole Fraction (Xco) for Variations in K^z and E-^^o 



H2O 



Eddy Diffusion Coefficient 



enrichment" 


4 X 10^ 




1 X 10^ 


4 X 10^ 




1 X 10^ 


0.51x 


1.16 X 10" 


10 


2.02 X 10-^° 


4.57 X 10" 


-10 


7.71 X 10-10 


Ix 


2.28 X 10- 


10 


3.96 X 10-10 


8.95 X 10- 


-10 


1.51 X 10-9 


2x 


4.57 X 10- 


10 


7.93 X 10-10 


1.79 X 10- 


-9 


3.03 X 10-9 


4x 


9.18 X 10- 


10 


1.59 X 10-9 


3.60 X 10- 


-9 


6.08 X 10-9 


8x 


1.85 X 10" 


-9 


3.22 X 10-9 


7.27 X 10- 


-9 


1.23 X 10-^ 


''Water enrichment factor [E-a^o) 


relative to 


a solar H2O/H2 


ratio of 9.61 x 10" 


4 


The Galileo entry probe 



measured a deep Jovian water abundance of 0.51x solar (IWong et alj|2004r). The observed mole fraction of 
CO in Jupiter's atmosphere is Xco = (1-0 ± 0.2) x 10-9 (|Bezard et al.ll2002h . 
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Table 3: Important Reaction Pathways for C0<=^CH4 in Jupiter's Atmosphere 



Reaction Rate Constant Reference 



R742 


H + CO 




HCO 


fco = 5.30 X lo-34e(-3™/^) 


Baulch et al. (19941 










= 1.96 X io-i3e(-i366/T) 


Arai et al. (1981) 


R743 


HCO 




CO + H 


fco = 4.18 X 10-9r-0-36e(-8294/T) 


reverse of R742 










= 1.54 X ioi2T-0-36e(-9290/T) 


reverse of R742 


R786 


H + H2CO 


— >■ 


HCO + H2 


fcR786 = 9.53 X io-i7Ti-90e(-i379/T) 


Irdam et al. ( 1993) 


R787 


Ha + HCO 


— >■ 


H2CO + H 


fcR,787 = 3.92 X 10-19T2-23e(-9082/T) 


reverse of R786 


R830 


CH3O 




H2CO + H 


fco = 1.40 X 10-6r-l-20e(-7800/T) 


Page et al. ( 19891 










fcoo = 1.50 X ioiiTe(-i28SO/T) 


Bezard et al. (20021 


R831 


H + H2CO 




CH3O 


fco = 1.28 X 10-33T-0.30e(3074/T) 


reverse of R830 










fc^ = 1.37 X 10-16rl-90e(-2006/T) 


reverse of R830 


R862 


H + CH3OH 


— i> 


CH3O + H2 


fcR862 = 1-135 X 10-22T3-«e(-3640/T) 


Jodkowski et al. (1999) 


R863 


H2 + CH3O 


-> 


CH3OH + H 


fcR863 = 1-77 X 10-22r3-09e(-3055/T) 


reverse of R862 


R858" 


H + CH3OH 


-> 


CH3 + H2O 


fcR858 = 3.321 X 10-10e(-2670/T) 


Hidaka et al. ( 19891 


R859 


H2O + CH3 


— >■ 


CH3OH + H 


fcR859 = 5.83 X 10-15Te(-15474/T) 


reverse of R858 


R150 


H + CH4 


— i> 


CH3 + H2 


fcRlSO = 2.20 X 10-20r3-00e(-4045/T) 


Baulch et al. (19921 


R151 


H2 + CH3 


-> 


CH4 + H 


fcRlSl = 8.62 X 10-24T3-57e(-2995/T) 


reverse of R150 



Note - Reaction numbers in bold refer to steps in the dominant CO— >'CH4 kinetic scheme. Rate constants for bimolecular 
reactions (fci for reaction i) and high- pressure limiting rate constants for termolecular reactions (fccxD) are in units of cm3 s~^. 
Low-pressure limiting rate constants for termolecular reactions (fco) are in units of cm^ s~^. For decomposition reactions (e.g., 
R743 and R830), fco is in units of cm3 s~^ and k^o is in units of s~^. Note that although the rate coefficients were internally 
reversed within our code, we provide empirical fits to the reaction rate coefficient expressions (valid for T = 400 to 2500 K) for 
reverse reactions as an aid to other investigators. "See note added in proof. 
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Table 4: Model results for variations in H + CH3OH reaction kinetics 



Eh 


2O 


RLS 


Rate Coefficients for H + CHaOIf Reaction 


Pathways" 


0.4- 


-3.4 


R863 


2.249 X 10-2lT3-2e(-"55/T) ^ 135 y. io-22rp3A^(-3640/T) 


9.41 X io-9e(-l2400/T) 


0.5- 


-4.8 


R863 


Jodkowski et al. figgg) Jodkowski at al. fl999) 

5.47 X 10-15Tl-24e(-2260/T) 2.28 X l0-i6rpl.24,^(-226Q/T) 


fcR860/A:R858 = 40 at 1000 K 
1.82 X io-8e(-l2400/T) 


0.6- 


-5.3 


R858 


Baulch et al. f2005), assuming A;R«fin /fcRSR9 = 24'' 
2.72 X io-i7r2e(-2273/T) g 33 X io-i7T2e(-2947/T) 


A:R86o/fcR858 = 40 at 1000 K 
1.70 X io-8e(-i2400/T) 


0.1- 


-1.2 


R858 


Li and Williams (1996), fcRaKn/AiRSK, = 0.43e674/T 
4.23 X 10-10e(-3516/T) i ig X 10-10e(-5549/T) 


fcR860/fcR858 = 40 at 1000 K 
5.08 X io~iOe(-i2400/T) 


0.6- 


-5.0 


R858 


Lendvav et al. ( 1997") Lendvav et al. C 1997) 
4.56 X 10-15Tl-24e(-2260/T) X ^Q-15rpl.2i^(-22m/T) 


Lendvav et al. (1997) 
1.52 X l0-8e(-l2400/T) 


0.7- 


-6.2 


R831 


Baulch et al. (2005). assumine k-RURn /knuKo = 4"^ 
2.82 X io-i7r2-ie(-2«o/T) y 04 X io-i8T2-ie(-2450/T) 


fcR86o/fcR858 = 40 at 1000 K 
2.94 X io-8e(-i2400/T) 



Tsang (1987), fcR86o/fcR862 = 4 fcR86o/fcR858 = 40 at 1000 K 

Note - "Rat e coefficients arc in units of cm'^ for bimolecular reactions. ''Based upon lCarvalho et al. I boosi ). " Based upon 
iTsand lll987l) . 
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latitude (deg) 



Fig. 1. — Eddy diffusion coefficient K^z at non-equatorial latitudes as a funct ion of latitude 
and pressure predicted using the theory of iFlasar and GieraschI (119771 . Il978l ) for thermally 
driven turbulent convection in a rapidly rotating system (see Eq. (jS])). 
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Fig. 2. — Carbon equilibrium chemistry for CH4/H2 = 2.37 x 10^^ ( jWong et al.l 120041 ) and 
H2O/H2 = 2.40 X 10"^ (2.5x s olar) in Jupiter's atniospher e, calculated using the NASA CEA 
code. Modeled after Fig. 17 in lFegley and Lodderd ( 1l994j ) for a slightly different atmospheric 
composition and revised thermodynamic parameters for H2CO and CH2. 
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uL 




X 

o 



10 ^^10 ^° 10 ^ 



10 ^ 10 ' 10 " 10" 
mole fraction abundance 



-7 



10^10" 0.01 



Fig. 3. — Vertical profiles for the CO, H2O, and CH4 mole fractions in Jupiter's atmosphere 
for an assumed water enrichment of Ix solar (wher e we assume the solar value is H2O/H2 = 



9.61 X 10 ^, iLodderj (120041 ): iLodders et al.l (120091 ) and for assumed constant K^z values of 



4 X 10^ (dotted line), 1 x 10^ (dot-dot-dot-dash line), 4 x 10^ (dash-dot line), and 1 x 10^ 
cm^ s~^ (solid line), as well as a Kzz profile determined by Eq. (dashed black line). The 
dashed gray line represents the CO abundance predicted by chemical equilibrium using the 
NASA CEA code. The circle with error bars represents the observed tropospheric CO mole 
fraction reported by lBezard et al.l ( 12002| ). Note that if the deep water enrichment is Ix solar, 
as in the models shown here, the CO observations are best reproduced for Kzz = 4 x 10^ 
cm^ s~^. Note also that water and methane are the dominant oxygen- and carbon-bearing 
gases, respectively, throughout Jupiter's deep troposphere. 
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Fig. 4. — Vertical profiles for CO and CH4 in Jupiter's atmospere for an assumed eddy 
diffusion coefficient K^z = 1 x 10^ cm^ s~^ (our nominal model) and for various assumed 
water enrichments (£'h2o) of 0.51 (GPMS value; dotted line), 1 (dot-dot-do t-dash line), 2 
(dash-dot line) and 4 times (solid line) the solar H2O/H2 ratio of 9.61 x 10""^ ( lLoddersll2004l : 
Lodders et al.ll2009l ). The corresponding gray lines show the CO abundance predicted by 



chemical equilibrium using the NASA CEA code. The circl e with error bars re presents the 
observed upper tropospheric CO mole fraction reported by iBezard et al.l (|2002[ ). Note that 
for our assumed nominal value of K^z = 1 x 10^ cm^ s~^, the CO observations are best 
reproduced for assumed global water enrichments of 2-4 times solar. 
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water enrichment factor 



Fig. 5. — The quenched CO mole fraction derived from our model (solid lines) as a function of 
the water enrichment (E'hjo) relative to the solar value of H2O/H2 = 9.61 x 10~^ for different 
values of (as labeled). The dashed vertic al line shows the water abundance measured 
by the Galileo entry probe (IWong et al.l 120041 ) . The shaded area i ndicates the observed CO 
mole fraction of (1.0 ±0.2) x 10~^ reported by lBezard et al.l (120021 ) and constrains the range 
of Kzz and E-^^q values in Jupiter's troposphere. The error bars for the K^z = 1 x 10^ 
cm^ s~^ solution represent the estimated errors due to uncertainties in the reaction kinetics. 
Higher Kzz values represent quenching and mixing from deeper in the atmosphere where 
CO is more abundant; higher water enrichments yield higher CO abundances at the quench 
level. 
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^Q-12 ^Q-6 

mole fraction abundance 

Fig. 6. — Vertical abundance profiles for oxidized carbon gases in Jupiter's atmosphere for 
Kzz = 1 X 10^ cm^ and H2O/H2 = 9.61 x 10~^ (Ix solar). The dotted gray lines indicate 
thcrmochcmical equilibrium; divergence from the equilibrium profiles (shown here for CO, 
CO2, H2CO and CH3OH) show where rapid vertical mixing and slow reaction kinetics drive 
each species to a constant quenched profile. For HCO, CH2OH, and CH3O, disequilibrium 
mixing occurs at much lower abundances (not shown). 
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Fig. 7. — Observed Jovian heavy element enrichments relative to solar (filled circles), as 
measured by the Galileo entry probe ( C as CH4, N as NH3, O as H2O, S as H2S, Ar, Kr, Xe; 
Wong et al.ll2004j : iMahaffy et al.ll2000l ) and infrared spectroscopy (P as PH3, see §2.31) . along 
with our kinetics-transport model results ( open circle) for 0_as H2O. Solar heavy element-to- 
hydrogen abundance ratios are taken from lLodders et al.l (l2009l ). The water enrichments are 
defined relative to a solar water abundance of H2O/H2 = 9.61 x 10~^ (see text for details). 



